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ABSTRACT: We propose that there exist generalized Seiberg-Witten equations in the Liouville 
conformal field theory, which allow the computation of correlation functions from the resolution of 
certain Ward identities. These identities involve a multivalued spin one chiral field, which is built 
from the stress-energy tensor. We solve the Ward identities perturbatively in an expansion around 
the heavy asymptotic limit, and check that the first two terms of the Liouville three-point function 
agree with the known result of Dorn, Otto, Zamolodchikov and Zamolodchikov. We argue that 
such calculations can be interpreted in terms of the geometry of non-commutative spectral curves. 
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1. Introduction 

Among two-dimensional conformal field theories, Liouville theory may be viewed as the simplest 
non-trivial theory with a continuous spectrum. Unlike minimal models, which exist for discrete 
values of the central charge c, Liouville theory can be defined for any real value c > 25, and ana- 
lytically continued to arbirtrary complex values of c. Like minimal models, Liouville theory is min- 
imal, in the sense that it is completely determined by symmetry and consistency conditions, once 
the spectrum of primary fields is given. This spectrum is specified by the set of allowed conformal 
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dimensions A > ^j- , with multiplicity one for each dimension. The Dorn-Otto-Zamolodchikov- 
Zamolodchikov formula for the three-point function is the unique solution of the crossing 

symmetry equations for this spectrum |3|,Q]. And with the conformal bootstrap method , it is in 
principle enough to know the three-point function for computing all correlation functions. 

The minimal nature of Liouville theory probably explains why this theory has many funda- 
mental applications, in particular in two-dimensional gravity and exactly solvable models of string 
theory. In these applications, the conformal transformations which are symmetries of Liouville 
theory are respectively interpreted as space-time diffeomorphisms and worldsheet reparametriza- 
tions. More recently, applications to four-dimensional gauge theories were initiated by the work 
of Alday, Gaiotto and Tachikawa [§J. The existence of the AGT relation between Liouville and 
gauge theories may seem surprising, because two-dimensional conformal transformations do not 
usually appear in gauge theories. But the gauge theories in question belong to a special class of 
supersymmetric theories, where such transformations do appear. 

Since Liouville theory is so fundamental, it is worthwhile to investigate the structure of the 
theory, and to look for alternative methods of computing its correlation functions. We are motivated 
by the following question: are the correlation functions of Liouville theory encoded in the geometry 
of some spectral curves? This is inspired in particular by the relations of Liouville theory with 
matrix models, whose observables are known to be encoded in certain algebraic curves. However, 
matrix model constructions cannot be directly relevant to Liouville theory, which has no integer 
parameter corresponding to the integer size of the matrix. Only particular correlation functions 
in Liouville theory, where a continuous parameter takes an integer value, are related to matrix 
integrals. (See Appendix for more details.) However, the Topological Recursion formalism [pi], 



which was originally developed in the context of matrix models in order to compute the observables 
from the geometry of the spectral curve, has a more general validity, and requires neither a matrix 
integral nor an integer parameter. 

In this article we will propose an alternative derivation of Liouville correlation functions, based 
on techniques inspired from Topological Recursion. The logic of the argument is: 

1 . The basic principle of the relation of Liouville theory with matrix models (and of the AGT re- 
lation) is that each correlation function of Liouville theory can be interpreted as the partition 
function of a specific model, 

z = mv ai (z i ,z i )\ , (i.i) 

so that the variables ai (the momenta) and Zj (the positions in the two-dimensional space) 
of our Liouville n-point function are interpreted as parameters of a model whose partition 
function is Z. 

2. Using the holomorphic stress-energy tensor T, we define a chiral (i.e. locally holomorphic) 
field J by 

T = - J 2 + QdJ , (1.2) 
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where Q is defined by c = 1 + 6Q 2 . We introduce the correlation functions 

W m (xi,---X m ) = JiXj)^ , (1.3) 

where we define the double brackets by 

P}=^(of[V ai (z l ,z l )^ . (1.4) 

Then the correlation functions W m obey Ward identities, which are analogous to the loop 
equations of matrix models. 

In order to derive the partition function Z from the study of the correlation functions W m , 
we need to assume the validity of a Seiberg-Witten equation of the form 

d 



da.,' 



F= / Wi(x)dx, (1.5) 



Bi 



where Bi is some contour of integration, and the free energy F is defined by 

F = logZ. (1.6) 

The Seiberg-Witten equation is inspired by matrix models and supersymmetric gauge theory, 
and it can be heuristically justified from the functional integral representation of Liouville 
correlation functions. We will treat it as a basic assumption, which will be justified if it leads 
to reasonable results. 

4. Finally, we can ask whether there exists a spectral curve C where the 1-cycle Bi, and the 
1-form W\(x)dx live. However, this "exact" spectral curve is not very useful or accessible. 
As we will perform calculations in a perturbative expansion of Liouville theory, we will 
rather discuss the "perturbative" curve which is associated to that particular expansion, and 
on which the leading term W-[ (x)dx of W\(x)dx lives. More precisely, we will consider 
the heavy asymptotic limit 

Q -> oo , r)i = — fixed , (1.7) 
and the corresponding perturbative expansion. This leads to the non-commutative curve 

2/ 2 + t (0) (x) = , [y,x] = l, (1.8) 
where t^°'(x) = lim Q 2 (T(x)} is a rational function of x with second-order poles at Z{. 

Q— s-oo 

Using the identification y = -S^, the equation of our curve can be viewed as a differential 
equation, namely the Fuchsian form of the Liouville equation with n sources. 

Our concrete aim is to rederive the DOZZ formula for the Liouville three-point function using 
these ideas. We will partly fulfill this aim, by rederiving the first two terms of the expansion of 
that formula around the heavy asymptotic limit. This will be done using the formalism of [^] 
for correlation functions associated to non-commutative spectral curves. (However, no previous 
knowledge of that article or of Topological Recursion is assumed.) 
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2. The spin one chiral field in Liouville theory 



In this section, we study the chiral field J defined by eq. (1.2). Some readers may object firstly 
that the solution of the quadratic equation ( jl.2| ) cannot be a unique, well-defined chiral field J, and 
secondly that if it was, it would be a current, so that we would be studying a free boson instead 
of Liouville theory. The answer to both objections is that J should be considered as a multivalued 
chiral field, which is locally holomorphic and has not only poles like the stress-tensor T, but also 
branch cuts. Due to these branch cuts, some of the Ward identities of the free boson theory fail, 
and J does not encode a symmetry. The relation between the stress-energy tensor T and the 
chiral field J is actually nothing but the operator version of the usual construction of the Virasoro 
algebra in terms of the Heisenberg algebra. 

2.1 Extended chiral algebra 

The symmetry algebra of Liouville theory is the Virasoro algebra of local conformal transforma- 
tions. This algebra is encoded in the operator product expansion (OPE) of the stress-energy tensor 
T with itself, 

T(x)T(x') = - § + 2 T(x>) + -^dT(x>) + 0(1) , (2.1) 

[x — x p (x — x y x — X 

where c is the central charge. The vertex operators V a (z, z) are defined by their OPE with the 
stress-energy tensor, 

T(x)V a (z, z) = ( Aa + -^—f] V a (z, z) + 0(1) , (2.2) 
\{x — z) z x — z oz ) 

where the conformal dimension A a associated to the momentum a is 

A a = a{Q - a) , (2.3) 

where we recall the relation between the central charge c and the background charge Q, 

c = 1 + 6Q 2 . (2.4) 

Now in the free boson theory the symmetry algebra is the affine algebra u(l), which is encoded in 
the OPE of the chiral field J with itself, 

_ l 

J{x)J{x') = ^—r + J 2 (x') + 0(x - x') . (2.5) 

(x — x'Y 

The enveloping algebra of this u(l) algebra contains a Virasoro algebra with central charge c, 



whose stress-energy tensor T is constructed from J by eq. ( |1.2| ). We have written the J J OPE up 
to its first regular term J 2 (x'), as this provides the definition for the J 2 term in eq. (L2). The field 
J is quasiprimary of conformal dimension one with respect to T, 

r < l)J < I '» = (^ + ^^ J M +0(1) ' (2 - 6) 
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The formal identity of our construction with the free boson theory shows that the set of OPEs (2. 1), 



d2.5|), ( |2.6[ ) defines an associative algebra of operators. 

The crucial difference between the Liouville and free boson theories appears in the OPE of J 
with a vertex operator V a (z, z), which in the free boson theory reads 



J{x)V a (z,z) 



-a 



(free boson) X — Z 



V a (z,z) + 0(1) . 



(2.7) 



This equation expresses the conservation of the momentum a which comes from the u(l) symme- 
try. This symmetry is absent in Liouville theory. Actually, vertex operators with momenta a and 
Q — a have the same conformal dimension A Q (2.3), so that they must be identified, as a conse- 
quence of the assumption that each representation of the Virasoro algebra appears with multiplicity 
one in the spectrum of Liouville theory: 



V a (z, z) = R(a)V Q - a (z, z) , 



(2.8) 



where R(a) is a reflection coefficient This reflection equation is clearly incompatible with the 
OPE ( |2"77| ) as it stands. Notice however that we can make sense of that OPE in the context of the 
idea of a multivalued field J. Looking for a J(x) V ai (zi, Zi) OPE which would be compatible with 
the known T(x)V ai (z{, z~i) OPE ( |2J| ) given the relation (L2) between T(x) and J(x), we find the 
two solutions 



j(x)V ai (zi,Zi) 



^9-V ai {z h -Zi) + Oil) 



(2.9) 



where we take ± to label sheets of the space where J lives. Such sheets are not expected to be 
globally defined, and the labels ± in this OPE make sense only in the neighbourhood of x = Z{. 

2.2 Seiberg-Witten equations 

We now look for Seiberg-Witten equations involving correlation functions of the field J. Since 



the field J is multivalued, its correlation functions W m (|L3|) and in particular W\ are multivalued. 

Zi we have 



Wi(x) 



And the JV operator product expansion (2.9) implies that in the neighbourhood of x 
+ 0(1), which leads to 



a * = -hl^ )dx - 



(2.10) 



This equation is of the type of a{ = f A _ W\{x)dx, which is usually associated to an equation of 



the type of eq. (1.5) provided the contour Ai is dual to the contour B{ of that equation. This 



suggests that eq. (|L5j) may hold, provided the contour Bi ends at z\. Then Bi must have another 
endpoint. The choice of this other endpoint is problematic, and we will circumvent the problem 
by considering 
quantity is 



d 
f)t\, 



j F instead of just jj^-F. The natural Seiberg-Witten equation for this 



d 
dan 



d 

da.,. 



F 



W x (x) + W x (x) ) dx , 



(2.11) 
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whose right-hand side is antisymmetric under i O j as it should be. 



It may seem that our choice of using the solution W\(^x ) and not W\(x ) is arbitrary - and 



incompatible with the reflection symmetry eq. (|2.8|). We will later justify this choice under some 
assumptions ( |Q| ) on the values of parameters. And the results which we will obtain using our 
Seiberg-Witten equation will turn out to be compatible with reflection symmetry. 



In Appendix BJJ , we will prove generalizations of the Seiberg-Witten equations, which involve 
higher correlation functions W m instead of F and W\. 



3. Ward identities 

The multivalued field J is of course expected to have multivalued correlation functions. However, 
this issue can be ignored for the purpose of deriving Ward identities. The issue will matter only 



when we start solving these identities in Section |3_3 



3.1 Non-perturbative Ward identities 



It is convenient to write identities not for the correlation functions W m (1.3), but for their connected 
components W m defined as 



Wi(x) 
W 2 {xi,x 2 ) 

W 3 (xi,X 2 ,X 3 ) 



W x {x) , 
W 2 (xi,x 2 ) 



+ W 1 (x l )W 1 {x 2 ) 



Oi - x 2 y 

W 3 (xi,X2,x 3 ) + ^Wi(xj) {w 2 (x k ,x e ) 

3=1 V 



(x k - Xif 



which are examples of the general definition 



(3.1) 
(3.2) 



+ \{Wi{x 3 ) 



(3.3) 



W m (I) 



e n 

pgP(/) Jep 



W {Jl (J)-S l 



\J\,2 



(Ax(j)y 



(3.4) 



where / is a finite set of variables of length |/|, P(I) is the set of partitions of / into nonempty 
subsets, and Ax(J) is the difference of the two variables belonging to a set J of size two. So, 
W 2 (x,y) has a second-order pole at x = y due to the J J OPE (2^), but this pole is explicitly 
subtracted in the above definition, so that W 2 (x, y) is regular at x = y. Like W m , W m is symmetric 
under permutations of its arguments (x\, ■ ■ ■ x m ). 

The idea which leads to the Ward identities is to express a correlation function of the type of 



in two different ways. On the one hand, the x-dependence can 



be fully determined by the axiom that T(x) is meromorphic with poles only at Xj and Z{, and that 
the behaviour of T(x) near these poles is determined by the TJ ( ^^ ) and TV ( |2.2[ ) OPEs. (In 
particular, T{x) is supposed to be regular at x = oo in the sense T{x) = o(l).) On the other 
hand, the T insertion can be written in terms of J using the definition (|L2|) of J (together with eq. 



(2.5) for the regularization of J ). 
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In the case m = 0, we thus want to express (T(x)}} in two different ways. On the one hand, 
meromorphy of T(x) leads to 



(T(x)j = t(x) , 



(3.5) 



where the double brackets are defined in eq. (1.4), and we define the function t(x) as 



A, 



OZi 



A2 



+ 



Pi 



X Z{ 



(3.6) 
(3.7) 



We will call the quantum accessory parameters. This is because the leading term of appears 



as an accessory parameter in the second-order differential equation ( p.33| ) (see for instance Q9p ), 
which is the heavy asymptotic limit of the first Ward identity. On the other hand, inserting the 



definition (1.2) of J into eq. (3.5) leads to the first Ward identity, 



Q-^W 1 {x)-W 1 (x) 2 -W 2 (x,x)=t(x) . 



(3.8) 



This is the m = case of the Ward identity, which is obtained an equality between two different 
expressions for yT(x) Ilj=i J( x j))' 



Q^-W m+1 (I,x) - lim W m+2 (I, x,x') + 2 —W m (I) 



t(x) +D X + J2 



a i 



(x - x'y 

m ^ 

W m (I) + Qj2 f v, Wm-iCT - {Xj}) , (3.9) 



(x — Xi) 

j=l J J J j=l ^ 31 

where we denote I = (x\, x 2 , ■ ■ ■ x m ), and we introduce the differential operator 

1 d 



i=i 



dzi 



(3.10) 



Rewriting these Ward identities in terms of the connected correlation functions W m defined by eq. 



(3.4) yields 







Q-^W m+ i(I, x) - W \ J|+i( J ' x)W m _\j\ +1 {I — J,x) — W m+2 (I, x, x) 



Jd 



6 m>0 t(x) + D x W m {I) + J3 



d w m (i)-w m (i -{ Xj },x) 



dxj 

3=1 J 



(3.11) 



whose special case m = coincides with eq. ( p.8[ ). Notice that for m ^ each term of the sum 
over the subsets J of I appears twice, as the J and I — J terms coincide. The Ward identities are 
special cases of the loop equations of [||] (see eq. (5.9) therein), so that we can use the formalism of 
that article when it comes to solving them. Our identities are more specific than the loop equations, 
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which involve undetermined polynomial functions. This will allow us to go further in the resolution 
of Ward identities and calculation of free energies than can be done in the general setting of . 
One may however wonder how this formalism can possibly help us compute the partition 



function Z. Using the Seiberg-Witten equations (1.5) we can deduce Z from W\, which can itself 



be obtained by solving eq. (3J3) if we know W2, whose calculation as the solution of eq. ( |3.11 ) 



with m = 1 requires the knowledge of W3, and so on. On top of this, each Ward identity comes 
with its own integration constant. However, it will turn out that in the perturbative expansion of 
Section [3.3[ , computing W m up to the order g involves only W\, W2, ■ ■ ■ W m+g . For example, eq. 



( p.8| ) will lead to a closed equation for the leading term of W\ , as the W2 term will be subleading. 
Moreover, the integration constants will be uniquely determined by simple analyticity conditions. 

In the case of the free boson theory, the field J(x) in W\{x) = (J(x)) is assumed to be 
meromorphic with poles only at Zj, and the behaviour of J(x) near these poles is determined by 
the JV (O) OPE. This leads to 



Wi(x) = (3.12) 

(free boson) X — ?*-■ 

i=l 

This can be generalized to ( J(x) n*Li J( x j) /> 



m 1 

W m+1 (I,x) = WL(ar)W- m (/)-y; ? ^-,W m ^{I - { Xj }) . (3.13) 

(free boson) j—\ ~ X 3' 

This implies the vanishing of all higher connected correlation functions, 

W m > 2 = 0. (3.14) 

(free boson) 

3.2 Global conformal symmetry, and the case n = 3 

Let us analyze how correlation functions W m behave under global conformal transformations 
z — > Hqrf • Symmetry under such transformations constrains how W m depends on (x\, ■ ■ ■ x m ) 
and {z\, ■ ■ ■ z n ), and if n < 3 these constraints can be interpreted as completely determining the 
dependence on (zi, ■ ■ ■ z n ). 

The consequences of global conformal symmetry can be neatly encoded in the behaviour of 
the chiral fields T(x) and J(x) as x — > 00, 

T{x) = °{v) ' J(x) = "f + ' (3 " 15) 

which follow from the TT and TJ OPEs, and the assumption that the stress-energy tensor T 
generates the conformal transformations. This asymptotic behaviour must also apply to correlation 
functions of T and J. For instance, in the free boson theory, the study of the asymptotic behaviour 
of W\(x) ( T3.12 ) as x — > 00 leads to the conservation of momentum, Y%=i ai = Q- 



The Ward identity for T CiS) implies t(x) = O (^) , which constrains the quantum accessory 



parameters /3j of eq. (3.6). In the case n = 3 the accessory parameters are completely determined 
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by such constraints, and we have 

3 



(n=3) (£(x) ^ X - Zi 



where we introduce the function 

3 

<p(x) = l[(x - Zi ) . (3.17) 



i=l 



More generally, consider the Ward identity ( j3-9[ ), which is a relation between two expressions 
for yT(x) Ilj=i J( x j)f- m particular, its right-hand side must behave as O (^) as x — >■ oo. 
Requiring the vanishing of the 0(^),0(^) and O(^) terms leads to three partial differential 
equations for W m . This leads to constraints on the dependence of W m on (x\, ■ ■ ■ x m ), whose 
general solutions are 

Wi(x) = — -)H(u), u=^^, (3.18) 

X — Z\ \X — Z\ X — Z 2 J [X — Zs)Z21 

W m>2 ( Xl , ■■■x m ) = f[ (— —) H(u u ■■■u m ), Uj = p- - Zl \ Z23 , (3.19) 

jJ[\ X j~ Z i Xj-Z 2 J [Xj - Z 3 )Z 2 l 

where H(u) and H(u\, ■ ■ ■ u m ) are arbitrary functions. These expressions are symmetric under 
permutations of (z\, z 2 , £3), although this is not manifest. 

In the case n = 3, equations ( |3.18 )-( 3.19 ) determine the dependence of W m not only on 
(xi, • • • x m ), but also on (z\, z 2 , Z3), since there are no coordinates Z{ apart from (z\, z 2 ,z%) for 
forming conformally-invariant cross-ratios ( Zl ~ Zl }* 23 . This allows us to express the action of the 
differential operator D x ( 3.10| ) on W m as 

= -S mtl ^ + f2^( f Xj) ~^ W m (x lr --x m )) . (3.20) 

(n=3) (p[X) ~[OXj \ [X — Xj)ip[x) J 



This can be used for eliminating the D x W m (I) term in the Ward identity (3.1 1), with the result 

d 

Q-^-W m+1 (I, x)~Y^ W\j\ + i(J, x)W m _\j\ +l (I — J,x) — W m+2 (I, x, x) 



dx 

Jcl 



S mQ t( x ) 6 ml Q f 9 ipjx^Wmjl) -y{x)W m {I -{ Xj },x) (321) 
(n=3) m ' m ' <p{x) ~[dxj (x — Xj)ip(x) 

3.3 Expansion of the Ward identities 

To directly solve the Ward identities would be difficult. Instead, we will use the perturbative expan- 
sion (1.7) around the heavy asymptotic limit. Actually, there are several possible expansions around 
this limit. Our choice of Q as the expansion parameter, and of r\i = ^ as the fixed quantities, will 
lead to simple expansions for the Ward identities, but not necessarily for their solutions. Another 
possible choice would be to take b such that Q = b + b^ 1 as the small expansion parameter, and to 
have ba>i fixed. These two choices are related by a perturbative resummation of the expansion. 
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We assume that the expansions of F = log Z, W m and t(x) take the following forms: 

oo 

F = Y^ Q 2 ~ 2g F^ , (3.22) 

9=0 

oo 

W m = Y,Q 2 ~ m - 29 W$ , (3.23) 

(3.24) 

The behaviour of W m corresponds to the assumption that the spin one chiral field J is of order 



9=0 

oo 

t(x) = Y J Q 2 ~ 29 t {9) (x) ■ 

9=0 



Q, which is consistent with its behaviour at infinity eq. ( |3.15| ) . This implies that W m (\L2 
is of order Q m . The behaviour of the connected correlation function W m is then suggested by 
its definition (ET4l) from W m . The behaviour of t(x) (|3T6|) is suggested by the behaviour of the 



(3.25) 



(3.26) 



conformal dimension A a (12.31), 



QH , 5 = 77(1-7/). 



At a given order g > 0, the Seiberg-Witten equation ( [2.1 1[ ) reduces to 

w[ 9 \i) + w[ 9 \x) ) dx 



_d_ 

drjj drj. 



) F Li 



It is straighforward to expand the Ward identity ( 3.11 ) in powers of Q, and the 0(Q 2 2g m ) term 
of that identity is the perturbative Ward identity 



JClh=0 



3=1 J 



This identity can be used for determining provided we akeady know W^,^ for any pair 

(</, mf) such that 

(</ < g and m' < m + g — g') or (</ = g and m' < m) . (3.28) 

Therefore, the Ward identities allow us to recursively determine all Wm . This is what we now 
illustrate at low orders. The first two orders (g = 0, 1) of the expansion of the first Ward identity 

(m = 0) are 



d_ 

dx 



d_ 

dx 



2W[°\x) ) Wl L, {x) - W^'(x,x) = t^(x) . 



,(i) 



(0), 



(3.29) 
(3.30) 
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The leading order (g = 0) of the second Ward identity (m = 1) reads 



' d_ 

dx 



2W[°\x) ) Wr(x,y) = D x Wr(y) + 



,(o) 



(o). 



8 w[ 0) {y)-W[" > {x) 
dy x — y 



(o), 



(3.31) 



Considering (x) and (x) as known, we can solve eq. for then eq. (^31]) for 



and then eq. ( |3.30| ) for W^' . We could in principle go on and use equations with higher m 
and g, so as to determine successively W^°\ W% , etc. 

Let us comment on the equation (3.29) for W[ . Introducing the function \£ such that 



,(i) 



(o) 



,(0) 



the equation for becomes a second-order linear differential equation for 

V + t (0) W = o. 



(3.32) 



(3.33) 



This differential equation illustrates the multivalued nature of the spin one chiral field J. In partic- 
ular, in the neighbourhood of x = Zi, we have tW(x) = + 0(— ), which leads to two 



solutions ^(x ) and *&('x ) corresponding to the two critical exponents rji and 1 — r\i, 



^(x) = {x-z i f i {l + 0{x-z i )) , y(x) = (x-z i ) 1 - TH (l + 0(x-z i )) 

AO) 



(3.34) 



The corresponding solutions for W± are the only solutions of eq. (3.29) which are meromorphic 
in the neighbourhood of x = zi, and they behave as 



+ 0(1) , W^Cx 



1 



x 



+ 0(1). 



(3.35) 



This closely parallels the "multivalued operator product expansion" J(x)V a (z, z) eq. ( |2.9| ). 

A generic solution wf^ of eq. (3.29), built from an arbitrary combination ^ = Y2± A^^aT), 
must have a branch cut ending at z\. We will use the requirement of meromorphy near x = z\ for 
defining particular solutions of Ward identities, thereby solving the problem of the initial condition 
in these first-order differential equations. Of course, the resulting solutions still depend on the 
choices of the point z,- t and of the critical critical exponent rji or 1 — rfc. These choices will be 
indicated by superscripts as in 



4. Perturbative calculation of the three-point function 

We now want to solve the Ward identities, in order to derive the Liouville three-point function, 
which in our notations is the n = 3 case of the partition function Z (IT). We will focus on this 
case because the Liouville three-point function is explicitly known, and given by a non-trivial yet 
reasonably simple formula. (See Section |A.2j .) 

In the n = 3 case, some simplifications occur in our Ward identities. First of all, we can use 
the form ( 3.21 ) of the Ward identities, where the D x W m (I) term has been eliminated thanks to 
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global conformal symmetry. Moreover, the form ( [3.16 ) for t(x), together with the behaviour ( 3.25 ) 
of the conformal dimensions A ai , imply that t(x) vanishes beyond the leading order: 



tlx) = Q 2 t (0) (x). 

(n=3) 



(4.1) 



In the rest of this Section, we will omit the indication "(n = 3)". 

We will compute the first two terms and of the perturbative expansion ( 3.22 ) of 
F = log Z. To compute higher order terms looks complicated although in principle doable. 

4.1 Calculation of F(°) 

The task of computing F^ ' now boils down to solving the second-order equation ( 3.33 ), deducing 
W[ 0) thanks to eq. ( p2| ), and using the Seiberg-Witten equation ( |3.26| ) for obtaining F^°\ 

The second-order equation ( 3.33 ) can be solved explicitly in the case n = 3, because it 
amounts to a hypergeometric equation. For convenience, we reproduce the equation: 



X Za 



(4.2) 



where we assume i ^ j ^ k, and we recall the definition 5i = r/,(l — rji). We have introduced a 
labelling (3.34) of particular solutions, based on their behaviour near the critical points zi. We will 
now assume that the parameters obey 



(4.3) 



so that the solution ty^x) is dominant with respect to the solution \P(ic). We will work with 
dominant solutions only, and simplify the notations by writing $(z) = *&(x). For example, the 
dominant solution near x = z\ is 



= (x - Z\) 



X — Z2 \ ( X — Z3 x ' ' 



Z12 



Zl3 



F ^12^123 - l;2r/i; 



(X - Zl)z23 

(x - z 3 )z 2 i 



where we use the notations 

viz = vi + V2-v3 , V123 = m + m + m ■ 

The behaviour of ty(x) near x = zj with j ' ^ i is 

*(x) = Cij (x - zjp (1 + 0((x - z 3 ) l ~ 2 ^ + Q(x - Zj))) 



(4.4) 



(4.5) 



(4.6) 



where the coefficient Cij is an element of the monodromy matrix of the hypergeometric equation, 
which can be deduced from eq. ( A..8), 



_ T(2 m )T(l - 2 V2 ) ^ V2 ( z 13 

Z23 



rfa? 3 )r(i-»&) 



(4.7) 
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This is all the information we will need on the solution \P(x). The corresponding solution for 



W- 



(0), 



i yx) will be denoted W±°^ (x) = — #: log \P(x). It is meromorphic in the neighbourhood of z% 



and has a branch cut on (zj, Zf.) where i / j ^ k. 



Using the perturbative Seiberg-Witten equation (3.26) at the order g = 0, we compute 



d 



dm) F 



log 



log 



-I.I 



(4.8) 



mm**) 

where we use the regularization which consists in ignoring the divergent factor (x 
ty(zi). (This regularization is justified in 



m 



d 

dm 



d 

dm 



.) Using the formula ( |4.7[ ) for dj, we obtain 

= 2(t, 2 - 771) log z 12 + 2(7,! +772-1) log ^ + log ^j 7 ^ 3 ? , (4.9) 



Z13 



7(2771)7(77^3) ' 



where the special function 7(2) is defined in Appendix A.l. This suggests that the free energy 



which we obtain is related to the known result of Dorn, Otto, Zamolodchikov and Zamolodchikov 



Fdozz eq- (EH b y 



? (o) 

DOZZ 



<0) 



+ F 



(0) 



(4.10) 



where the bar indicates the replacement Zi — > z\. That our ft ) depends only on zi and not Z{ is 
an inevitable consequence of using locally holomorphic quantities W m , which is itself necessary 
for the integral in the Seiberg-Witten equation ( 2.1 1| ) to be independent from the precise path of 



integration. It is actually possible to derive the full dependence of Fqq ZZ on z\ and z%, by con 



sidering single- valued solutions of eq. ( 3.33 ) instead of our locally holomorphic solutions Vl/(x) 
[g]. However, that approach has been used only in the heavy asymptotic limit, which is the leading 
order of our expansion, and it seems difficult to generalize it to higher orders. 



4.2 Calculation of and the Gelfand-Dikii operator 



We will now solve equation ( |3.31[ ) for W 2 - This depends on the choice of a solution for W^' . We 



(o) 



will take an arbitrary solution = —%■, built from an arbitrary solution \l/(x) of the second- 



order equation (J4.2[). Let us introduce differential operators 



V x = d x + 2 



ty'(x) 
*(x) 



Vt = d x - 2 



*'(x) 
*(x) 



Using these operators, and the simplification ( J3.2QD which occur in the case n 
identity ( [3.31 ) can be rewritten as 

<p(v) 



v., 



(x - yf 



2\v!?\x,y) 



y)(f{x) 



(4.11) 
3, the Ward 

(4.12) 



One may wonder whether V x V y W2°\x, y) is really symmetric under x <-> y if we compute it 
using this equation. Checking this involves the third-order differential operator 



G x = V x d x Vt = dl + At^(x)d x + 2(t(°))'(x) . 



(4.13) 
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This is called the Gelfand-Dikii operator (see for instance Jl0| ] in [11]), and it has the property 

¥>0) „ <p(v) 



G, 



_ Q 

(y-x)(p(y) y (x-y)ip{x) 



(4.14) 



which follows from the fact that ip 2 t^ is a polynomial of degree two. This property directly leads 
to the expected symmetry of V x S/ y W^ ) \x, y). Notice that the Gelfand-Dikii operator appears 



in the classical higher equation of motion of Liouville theory [12], which in our notations reads 
G x ^{x) 2 = 0. More generally, the classical higher equations of motion of order m + 1 reads 

m 

ILL 



(d — 2j-^) ) ^ m = 0, and we expect that the corresponding differential operator appears 



in the calculation of W, 



(o) 



The solutions of eq. (4.12) are 



W^\x,y) 



2(x-y) 2± 2V 2 (x) v yj a i-y 



v(0 ' 



(4.15) 



(0) 



where a is an a priori y-dependent integration bound. This bound is constrained by requiring W 2 
to be a solution of eq. ( |4.12 ) with x -<->• y. This can be shown to hold if a = Zi for some i, provided 
$l{zi) = 0, which itself holds if ^Riji > 0. 

It remains to compute W% (x, x), which appears in the equation ( 3.30 ) for This is not 

completely trivial, because the regularity of W^\x,y) at x = y is not manifest in our formula 
( FBI ). We obtain 



wj,°\x,x) = ^V x 



Or. 



a i-X \^(x)ip{i) 



x — a 



where we define i^\x) from t^\x) in eq. (4.2) by 

(0) 



+ V x d x \og^{x) + 2t^(x) , 
(4.16) 



(4.17) 



To conclude this section on W^' , notice that (x, y) 



i x ^ y yi coincides with the Bergman 
kernel of [Q] (Section 4.3 therein), where many of our equations can be found in a more general 
form. 



4.3 Calculation of 



Let us proceed to the calculation of the solution w[ (x) of eq. ( |3.29[ ) which is holomorphic 



near x = z%, just like W^°\x) is meromorphic near x = Z{. This solution is formally written 



as W\ (x) = V~ W\ (x,x), where W 2 (x,x) is obtained by using the solution ty(x) and the 



r(0). 



AO) 



integration bound a = Z\ in eq. ( |4.16| ), and the operator Vj is obtained from V x ( [4.11D by using 
that same solution ^(x). We thus find 



W{ 1] {x) = i-d x I(x) 



1 



1 



2(x 



+ d x \o g y(x) + 2Vt 1 i ( -°\x) 



(4.18) 
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where we define 



(4.19) 



Each one of the last three terms of W± (x) is meromorphic near x = z% with a first-order pole at 
x = Zi (this requirement fixes the ambiguity in V~H^ (x)), and the sum of the residues vanishes. 

The first term \d x I(x) is holomorphic near x = z%. 

We now use the perturbative Seiberg-Witten equation ( 3.26| ) at the order g = 1, and compute 
. This involves integrating (x). The only term in W± (x) which does not have an obvious 

primitive is the last term, namely V^ 1 ^ (x). However, it turns out that it is easy to write a 

primitive for the sum VrH(°\x) + VjH^ix) which appears in w[ (x) + W^{x), 




where the constant Wij is defined as the Wronskian 



(4.20) 



(4.21) 



This leads to 



1 d _ ) F (i) 

dm dr]i 



I , , 1 . ^(x) 2 2*(x)*(s) 



2^) + 2^,_- 



+ 



(* i) • 



(4.22) 



We will compute this expression by evaluating each term. While the total expression is finite, 
individual terms can be divergent. While the total expression cannot depend on {z\ , Z2,z%) due to 
global conformal symmetry, individual terms can. 

To compute I{zk) with k = i or k = j, notice that the factor ip(x) in the integral expression 

i i 

( |4.19[ ) for I(x) obeys (p(zk) = 0. This allows us to replace and Vl/(£) factors with their leading 
behaviours near £ = z^. Using the integral formula (|A.6|) we obtain 



I{ Zi ) = i>(2 m )-i>{l) , I(x) 



log 



+ i/>{l-2r) j )-il>{l) + 0{x-z j ), (4.23) 



where the special function -0 is defined in eq. (4.3). The contribution of the term | log — fr~ i n e 1- 



( J4.22 ) is easily computed using the behaviours (334) and ( |4.6[ ) of ^f(x) near x = and x = zj 
respectively. Let us now study the last term in eq. ( |4.22 ). The integration constant in V^ 1 i^(x) 
is determined by the requirement that this function be meromorphic near x = Zi, and this leads to 



v;H^(x) 



fy(xf 



(4.24) 
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where the power-like divergence at £ = Zi is regularized by analytic continuation in the parameter 



n, 



. We thus have H|*Mvri£(o)(x) 



0. On the other hand, V i 1 i^°\x) has a logarithmic 



divergence at x = zj, which cancels the logarithmic divergences of the other terms of eq. (4.22), 
and which we will regularize by introducing a factor of (£ — Zj) e i in the integrand. Keeping this 
regularization implicit, we have 



*(aj)*(x) 



W; 



vr^(x) 



/* J di * (£) H (°) (£) definition 



E(r]f 2 ,r] 12 3 — l, 277i 



(i,j,fc)=(l,2,3) 

(«l,«2,Z3)=(0,l,Oo) 



(4.25) 



So, by sending (zi, 22, 03) to (0, 1, 00), we defined a function -E(a, 6, c) which depends only on 
(i]iiV2, %)> although it is more convenient to take (r]f 2 , ??i23 — 1, 2?7i) as its arguments. 



We know the coefficient Cij (4.7), and we can compute the Wronskian Wij (4.21) with the help 
of the transformation property (|A.8j) of the hypergeometric function, 



W!2 



r(2r/i)r(2r 72 



-z 



m +772—1 
12 



£13 
Z23 



771-772 



(4.26) 



This leads to the integral representation ( |A.21 ) for the function E(a, b, c). Appendix A3 is devoted 



to computing that integral. Combining the result with the other terms of eq. (4.22), we obtain 



1 °-) F m 

di] 2 drji 



log 



7(27/1)7(77^ 



23/ 



+ (vh - 5) tyfais) + ^(1 - vh)) - (vk - h) tyfaaa + ^(1 - 4*)) 



+ (2t? 2 - i) (V(2// 2 ) + ^(1 - 2772)) - (27?! 



(^(2 ?7 i)+V(l-2 ? 7i)) . (4.27) 



This suggests that our is related to the known result F^ ozz eq. ( A.19| ) by 

_ 77(1) j_ p(l) _ 2F (1) _ 



7?u; _ 77(1) 1 77(1) 



(4.28) 



5. Conclusions 

5.1 The results and their possible generalizations 

We have shown that our results for the first two terms F^ ' and F^ of the (logarithm of the) Li- 
ouville three-point function are consistent with already known results. One may object that our 
Seiberg-Witten equation ( 3.26| ) did not allow us to compute F^ 9 \ but rather differences of deriva- 
tives ( -T^p — F^ 9 \ thereby leaving a large indeterminacy. However, this indeterminacy is 



lifted if we assume that the reflection symmetry (2.8) holds, even if we do not know the reflec- 
tion coefficient, as explained in Appendix |A.4| . Therefore, we can compute and F^ up to 
unimportant normalization constants. 

We expect that all F^ 2 ) can similarly be computed, and agree with the known results from 
Dorn, Otto, Zamolodchikov and Zamolodchikov as in the cases of and F^\ eq. ( gig ) and 
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eq. (4.28). Thus, we expect the non-perturbative relations 



Fdozz = F + F or equivalently Z DO zz = \Z\ 2 . (5.1) 
This is analogous to the holomorphic factorization of black hole partition functions into topological 



string partition functions found by Ooguri, Strominger and Vafa [13], in the sense that our calcula- 
tion of Z involves integrals of locally holomorphic quantities - that is, topological objects in two 
dimensions. 

Completely solving Liouville theory would mean computing not only the three-point function, 
but also general n-point functions, which do not holomorphically factorize. Such n-point functions 
can however be reduced to locally holomorphic quantities, by using their decomposition into com- 
binations of three-point functions and conformal blocks. In order to compute an n-point conformal 
block, we would need to deal with the existence of n quantum accessory parameters ( |3.7[ ), of which 
n — 3 are not fixed by conformal symmetry. However, an n-point conformal block depends not only 
on the parameters of the n involved fields V ai (zi), but also on n — 3 extra parameters, namely the 
conformal dimensions in the internal channels. These extra parameters are expected to be related 
to the quantum accessory parameters. 

We expect our methods to work not only in Liouville theory, but also in other conformal field 
theories. We certainly need these theories to be non-rational, so that we can write Seiberg-Witten 



equations like eq. (|L5|) using continuous parameters like a; L . It seems possible to generalize our 
construction to conformal Toda theory, although this will not necessarily lead to the explicit com- 
putation of the three-point function, which is a difficult problem already at the leading order [jjjp. 
It would be interesting to generalize our construction to Wess-Zumino-Novikov-Witten models. In 
such models there are meromorphic spin one currents J a {z), and one may be tempted to use them 
in Seiberg-Witten equations as analogs of our chiral field J{z). But a preliminary analysis shows 
that we actually need multivalued spin one chiral fields, which do not coincide with J a (z). 

5.2 The non-commutative spectral curve 

Here we discuss the interpretation of our calculations in terms of the spectral curve whose equation 



is eq. (1.8). This interpretation is unfortunately not precise enough for justifying the form of the 



Seiberg-Witten equation ( |2.11| ). A comprehensive geometrical framework for non-commutative 
spectral curves has been developed in but it makes important technical assumptions which our 
curves do not obey. 

The spectral curve is supposed to be determined not only by the theory under consideration, 
but also by the particular expansion we want to use for the correlation functions. The leading 
order of the expansion is supposed to provide the equation of the curve, and the higher orders 
are supposed to be computed from certain geometrical quantities of that same curve. In principle, 
we may refrain from doing any expansion and consider the associated "non-perturbative spectral 
curve", but determining this curve would be equivalent to computing the partition function, and the 
curve would therefore be useless as a computational tool. So we want to consider an expansion of 
the partition function in terms of one or several parameters. The more parameters, the simpler the 
curve, and the more complicated the expansion. 

In the case of Liouville theory, there is only one natural expansion parameter - the central 
charge, and we have demonstrated that the corresponding curve and the associated quantities wi S ' 
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are tractable. We could nevertheless perform a further expansion with respect to a second parame- 
ter. In particular we could introduce a parameter h such that we would have [y, x] = h instead of 



[y, x] = 1 in the equation ( L8) of the spectral curve. The resulting spectral curve, which would by 
definition be obtained in the limit h — >■ 0, would then be commutative. But the expansion would 
then be a double series in Q 2 and h, which is an unnecessary complication as far as the calculation 
of the free energy is concerned. 

So, we consider the non-commutative spectral curve which is defined by the second-order 
differential equation ( |1.8| ). By analogy with the corresponding commutative spectral curve, the 
existence of two linearly independent solutions of that equation can be interpreted as the existence 
of two sheets above a point x. Since the solutions can be linearly combined, the definition of 
such sheets is rather arbitrary. However, we have found that near each singularity x = Z{ of 
the curve, it is possible to define two priviledged sheets corresponding to the two solutions with 
diagonal monodromy transformations, so that the corresponding correlation functions W^°\ l x) 
are meromorphic near x = Zi. Such priviledged solutions play a crucial role in lifting ambiguities 
in Seiberg-Witten equations and in solutions of the Ward identities. But their existence relies on 
the particular form of the spectral curve, which is why our calculations are difficult to interpret in 
terms of a general theory of non-commutative spectral curves. 

Such an interpretation would nevertheless be very interesting, because we might then access 
such a powerful result as the invariance of the free energy F under symplectic transformations of 
the spectral curve, which is known to hold in the case of commutative spectral curves [0]. This 
could then lead to simple proofs of relations between apparently different theories, such as the 
AGT relation. 

A. Useful technical results 
A.l Special functions 

We assume that Euler's Gamma function T(x) such that F(x + 1) = xT(x) is known. This function 
appears in the integral 

o V ' T(a + (3 + 2) 

The function j(x) is defined as a ratio of two Gamma functions: 

7(*) = Jr^K ■ (A.2) 



The function tp{x) is defined as 



T(l-x) 



i>{x) = ^logr(z) . (A.3) 



This implies in particular 



d 

— log 7(2) = ip(x) +ip(l -x) . (A.4) 
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The function ^(a;) has the series representation 



E 

n=0 



(x + n)(y + n) 



x-y 



(A.5) 



The function ip(x) also has integral representations, 



du- 



1 



u-l 



^(A + 1)-V(l) 



du- 



ll 



u-l 



loge + ^(l)-V(-A) + o(l). (A.6) 



We review some well-known properties of the hypergeometric function. First of all, we have the 
series representation 



F(a, b; c; x) = 

We have the transformation property 

T(c)T(c-a-b) 



r(c) 



r(o)r(6) 



E 

ra=0 



T(a + n)T(b + n) n 
n!T(c + n) 



(A.7) 



F(a, b; c; x) 



+ (l-x) 
The last two formulas imply 

oo 

E 



r(c-a)r(c-6) 

c -a-b r ( c ) r ( a + ft - c ) 



F(a, b; a + 6 — c + 1; 1 — x) 

F(c — a, c — b;c — a — 6 + 1; 1 — x) 



r(a)r(6) 

r(a + m)r(6 + m) r(a)r(6)r(c - a - b) 



(A.8) 



(A.9) 



m=0 



(A. 10) 



mW(c + m) T(c — a)T(c — b) 

Another useful transformation of the hypergeometric function is 

F(a, b; c; x) = (1 - x) c_a - 6 F(c - a, c - 6; c; x) . 
A.2 Expansion of the Liouville three-point function 

We recall the Dorn-Otto-Zamolodchikov-Zamolodchikov formula [|l|, |2|] for the three-point func 
tion Zdozz = (V ai (zi)V a2 (z 2 ) V a3 (2:3)) of Liouville theory at central charge c: 

ZDOZZ = \Z12 | 2{A "3 " A °l ~ AQ2 } I Z 23 | 2{Aq 1 _ Aq 2 - Aq 3 > I Z13 | 2(A °2 - Aq 1 - A °3 ) 

xAm ELi« T fe (2a 1 )T 6 (2a 2 )T b (2a 3 ) 



(A. 11) 



T b (ai23 - Q)^&(ai 2 ) T fe(«23) T 6( a 13) 

where c = 1 + 6Q 2 and Q = b + b^ 1 . Here A and are Oj-independent normalization constants, 
which can be determined neither in our approach (see Appendix A4 ) nor in the conformal bootstrap 
approach. The special function T&(x), which depends on the parameter b and is invariant under 
b — > b~ l , is defined by 



log T fe (x) 



00 dt 
T 



Q 

X 

2 



sinh 2 [ S - x ) £ 



sinh y s mn ^ 



(A. 12) 
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Let us compute the expansion of Fdozz = log Znozz around the heavy asymptotic limit (1.7). 
To begin with, we compute the expansion of logT b(Qv) w i tn V fixed. In this computation, we 
neglect the ^-independent terms as well as the terms in (rj — |) 2 , because such terms contribute 
to the normalization constants A, fi in Zdozz- Up to such terms, we find that log T{,(Qr/) has an 
expansion in powers of Q~ 2 , where the leading term is of order Q 2 , as in eq. (3.22), 



\ogT b (Q V ) = Y,Q 2 ' 2g n {9) {v) 

g=0 



(A. 13) 



We find 



u^°\r]) = v{rj) 



u {2) ( V ) 



-^)-(2l(l-l)-sKW 



where we use the function 



v(r)) 



log j(t)dt. 



(A. 14) 
(A. 15) 
(A. 16) 



(A. 17) 



In our expansion, the Zj-dependent prefactors of Zdozz contribute only at the leading order, and 
we have 



7 (o) 

DOZZ 
i(<?>l) 



5f 2 log \z 12 | 2 - (523 log \ z 23 I 2 - $13 l°g 1*13 1 2 , 



A DOZZ 

where we define the functional 

T {m} (f) = f(2 m ) + f(2 m ) + f(2 m ) - f(rj* 2 ) - f(4 3 ) - f( V 2 3 ) - f( m23 - 1) 

A.3 Calculation of integrals involving hypergeometric functions 

The following integral appears in the calculation of more precisely in eq. ( 4.25 ): 

r(a)r(6)r( c - a)r(c - 6) 



(A.18) 
(A. 19) 

(A.20) 



E(a, b, c) 



4r(c) 2 r(c -a- b)T{a + b - c + 1) 



d^e i+c - x o--o 



£2+a+b—c 



(c- l) 2 (a + b-c) 



(a - bf 



(A.21) 



where we introduce two regulators ei and e 2 . Let us decompose it into three terms, as suggested by 
the last factor: 



E(a, b, c) = -- [(c - l)E l + (a + b - c)E 2 + (a - b)E 3 ] 



(A.22) 



The regulators e\ and e 2 are needed for the terms E\ and E 2 respectively, and will be set to zero 
in the other terms. In each of the three terms E±, E 2 and £"3 we expand one factor F(a, b; c; £) in 
powers of £, and the other factor F(a, b; c; £) in powers of (1 — £), using the formulas ( |A.7| ) and 
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( |A.8| ) for the hypergeometric function. Then we perform the integral over £ using eq. ( |A.1| ). The 
result is 



^ Y{c-l + m + ei) (T{a + m)T(b + m) T{a + n)T(b + n) 



1 ^ 7 m\T(c + m) \ T(a)T(b) n\T(a + b + m + n + ei 

m,n=0 

T(c — a + m)T(c — 6 + fn) T(c — a + n)T(c — b + n) 



T(c — a)T(c — 6) n!T(2c — a — 6 + m + n + ei) 



(A.23) 



where in the second term we used the transformation QA.10| ) before expanding the hypergeometric 
function in powers of £, 



a + b — c ^\ ( T(a + m)T(b + m) T(a + n)T(b + n) T(a + b — c + n + €2) 
2 = WnYP(h\ 2^ 



T(a)T(b) ^ \m\T(a + b + m + n + €2) n\ T(a + b- c+l + n) 

' m,n=0 v v ' x ' 

T(a + m)T(b + m) T(c - a + n)T(c - b + n)T(n + e 2 )\ (A 24) 
m!T(c + m + n + €2) T(c — a — 6 + 1 + n) n! / 



fc — a ^> r(a + m)r(6 + m) 
3 = T(n\T(h\ 2-~> 



T(a)T(b) ^ ml 

K ' v ' m,n=0 



r(a + n)T(6 + n) T(c - a + n)T(c - b + n) 



n\T(a + 6+ l + m + n) T(c + 1 + m + n)r(c — a — b + 1 + n) 



(A.25) 



These sums can be performed, and expressed in terms of the ip function eq. (A3), 



Ex = ip(c - a) + ip(c - b) - V>(a) - ^(6) , (A.26) 
2 

E 2 = — + 2^(a + b - c) + 2^(c - a - 6) - -0(c - a) - Mc - b) - Ma) - Mb) + 0(e 2 ) , 
£2 

(A.27) 

£ 3 = ^( a ) + ^(c - a) - y>(6) - V(c - b) . (A.28) 
We now give some details on the computations which lead to these results. 

Computation of E3. In each term of £3, we use the formula ( |A.9| ) for performing the sum over 
m. The sum over n is then of the type ( \S), and directly produces the tp functions. 

Computation of E%. In each term of E2, we use the formula ( |A.9| ) for performing the sum over 
m. Then we separately consider the term n = and the rest of the sum. The term n = diverges 
as €2 — > due to a T(e2) factor, and we have to expand the remaining Gamma functions in powers 
of €2 up to 0(e2), using ^-fejiy^ = 1 + etp(x) + 0(e 2 ). The n > 1 terms all have a finite limit as 
62 — > 0, and their sum is of the type (|A.5[). 



Computation of E\. In both terms of E\, we use the formula ( |A.9| ) for performing the sum over 
n. The rest of the calculation is analogous to the calculation of E2, with the m = term diverging 
as ei — y 0, and the rest of the sum being of the type ( |A.5| ). Combining the two terms of E\, the 
divergent O(-) terms cancel. 
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A.4 Reflection symmetry 

Let us show that the knowledge of the derivatives 

A, - (4- 







dcti daj 



(A.29) 



which appear in our Seiberg-Witten equations (|2.11| ), together with the assumption of reflection 
symmetry (2J5), determine F up to normalization constants. At the level of the free energy, the 
reflection symmetry implies 



F(ai,a 2 , ■■■) = F(Q - a x , a 2 , ■ ■ ■ ) + r(ai) , 
for some unknown function r(a). This implies in particular 

fi 2 {Q - ai,a 2 ,- ■ ■ ) + fi 2 (at,Q - a 2 , ■ ■ • ) = r'(ai) - r'(a 2 ) 



(A.30) 



(A.31) 



This equation allows us to check that f\ 2 is compatible with reflection symmetry, and to determine 
r(a) up to a linear function of a. We can then compute -^F thanks to 



2- — F = f\ 2 (a\,a 2 , ■■■) - f\ 2 (Q - a%, a 2 , ■ ■ ■ ) + r'(ai) 

OOL\ 



(A.32) 



This determines F up to a linear function of <x\. By permutation symmetry, the ambiguity in F 
amounts to the choice of two constants, which correspond to the parameters A and [i in eq. ( |A.ll[ ). 



B. Peripheral topics 



B.l Higher Seiberg-Witten equations 

Here we argue that the following higher Seiberg-Witten equation holds for any m > 1: 











da 4 dan 



W m {I) 



W m+1 (I,x) - S m> i 



(x — X\) 2 



dx 



(B.l) 



where / = (x\, ■ ■ ■ x m ). These equations are generalizations of the Seiberg-Witten equation for 
F (2.11), if we formally define Wq by Wq = F. In this sense, the existence of higher Seiberg- 
Witten equations, which we will prove in the case n = 3, lends support to the validity of the 
Seiberg-Witten equation for F. 

There is however a substantial difference between the Seiberg-Witten equations for F and 
for W m . The former equation involves a specific choice of a combination of solutions W\{x ) + 

3 + 

Wi( x ), in other words a specific determination of the integration variable x. The latter equations 



are however insensitive to this choice, as is indicated in eq. (B.l ) by the absence of superscripts for 
the variable x. 



To prove the higher Seiberg-Witten equations ( JB. it is enough to show that they are com- 
patible with the Ward identities ( 3.11| ). This is because the correlation functions W m are unique 



solutions of these identities, as shown perturbatively in Section p.3| . The compatibility is easy to 
prove (not forgetting that D x in eq. ( 3.1 1| ) acts on the integration bounds in eq. ( |B. lh ), except in 
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the m = Ward identity, which involves the function t(x). Let us apply the differential operator 
Z5q" — to the m = Ward identity eq. (^8), assuming the higher Seiberg-Witten equations 
hold. The integrand in the resulting equation can be simplified with the help of the m = 1 Ward 
identity, and we are left with the equation 

° *(*) = T-^-vI ~ T-^-vI + 2 f J (D x W x {y) + ±-L- Wl (y)) dy . 



datj doLi) (x - zj) 2 (x - Zi) 2 J z . \ dy x - y 

(B.2) 

To prove this equation for n > 4 would require some information on the quantum accessory 
parameters which appear in t(x) (^^), and this is outside the scope of this article. However, in 
the case n = 3, the accessory parameters are known. In this case, using the explicit expression 
( 3.16 ) for t(x) and the formula ( 3.2C| ) for D x W\{y), the above equation is found to hold, provided 

Vfc, Wi(y) = -^^ + o(^—) . (B.3) 

y~z k y - z k \y - z k J 

This condition is obeyed by almost all solutions W\{y), provided K(afc — S) < 0. This is easily 

understood in the heavy asymptotic limit, where our condition on a k becomes 3?% < \. Then 

all solutions ^f(y) of the linear differential equation ( |3.33| ) behave as ^f(y) = 0((y — Zk) r,k ), 

y~z k 

k— 

except the particular solution y ) which corresponds to the critical exponent 1 — 77^. Thus we 
have Wi (y) = — ZS y \ = + o ( — - — ) for almost all solutions. 

Our higher Seiberg-Witten equations are special cases of equations which describe the be- 
haviour of W m under arbitrary deformations of the spectral curve, and which can be found in [||] 
(equation (6.48) therein). 

B.2 Comparison with matrix models 

An n-point function (ni=i ^ -^)) °f Liouville theory at central charge c = 1 + 6Q 2 = 
1 + 6(6 + can be expressed as an iV-dimensional integral provided [Q] 

n 

^ai = Q-Nb, (B.4) 
i=l 

where N is a natural integer. In this case, the functional integral description of Liouville theory 
implies that (JlILi Va t {zi, z~i)) has a pole, whose residue is (up to unimportant factors) 

„ N n 

z N =n m" 4qiqj j dY n n - z^- 1 ^ . ^ 

i<j " a=l i=l 

where we introduce the following measure on the iV-uple of complex variables Y = (yi, • • • yjv") 

N 

dY = H\y ab \- 4b - 2 Hd 2 ya . (B.6) 

a<b a=l 

We can interpret {y a } as the eigenvalues of a complex matrix M of size N, and dY as a measure 
on the space of such matrices (the "beta ensemble"). This is the basis of the proof of the AGT 
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relation by Dijkgraaf and Vafa in this case Q15|]. The free energy Fjf = log Zn is moreover known 



to obey Seiberg-Witten equations. This was proved perturbatively in [16], elaborated upon in [|1 
and checked in the related case of the Gaussian potential in [|T|]. 

It is possible to apply our formalism in this case, and the spin one chiral field is then 

1 - 1 

J(x) = Tr — = V . (B.7) 

y ' x-M j^x-ya 

The hypergeometric solution ^f(x) (|4.4[) of our second-order equation reduces to a Jacobi polyno- 



N 



mial of degree N, 

3 

= JJ(x - JJ (x - Sa ) , (B.8) 

i=l a=l 

whose zeroes s a must obey 

E^^ + E^^ = °- ( B - 9 ) 

Actually, the solution Vl/(x) has a diagonal monodromy around not only Z\, but also Z2 and Z3, and 

123 

we have \l/(xl = ^(x) = \l/(x). It is natural to choose this solution in computing the correlation 
functions Wm \ which will then be meromorphic functions. On the other hand, the subdominant 

1— 2— 3— 

solutions x ), \&( x ) and ^f(x) remain different from each other, and their monodromy prop- 
erties do not simplify. 

So, in the special case (BA), the spin one chiral field J is meromorphic, as in the free boson 
theory (which is the case N = 0). It may seem that having polynomial solutions of the second- 
order equation brings important simplifications. And indeed, relatively simple formulas for the 
correlation functions Wm can be obtained. However, these formulas depend not only on the 
physical variables Zi, but also on the auxiliary variables s a . And it is not easy to eliminate the 
auxiliary variables in the resulting expression for the free energy F^K Therefore, when it comes 
to computing the free energy, the special case is not necessarily easier than the general case. 



Acknowledgments 

We are grateful to Michel Bergere for collaboration on related topics. We thank Ivan Kostov, Vin- 
cent Pasquier and especially Michel Bergere for helpful comments on the manuscript. L. Chekhov 
is grateful to the Russian Foundation for Basic Research for support (Grants Nos. 11-01-00440-a, 
ll-01-12037-ofi-m-2011, and 11-02-90453-Ukr-f-a), to the Grant for Supporting Leading Scien- 
tific Schools NSh-46 12.2012. 1, and to the Program Mathematical Methods of Nonlinear Dynam- 
ics. 



References 

[1] H. Dorn, H. J. Otto, Two and three point functions in Liouville theory, Nucl. Phys. B429 pp. 375-388 
(1994), hep-th/9403141 



-24- 



[2] A. B. Zamolodchikov, A. B. Zamolodchikov, Structure constants and conformal bootstrap in 
Liouville field theory, Nucl. Phys. B477 pp. 577-605 (1996), hep-th/ 95 613 6 

[3] J. Teschner, On the Liouville three point function, Phys. Lett. B363 pp. 65-70 (1995), 

hep-th/9507109 

[4] J. Teschner, A lecture on the Liouville vertex operators, Int. J. Mod. Phys. A19S2 pp. 436-458 (2004), 

hep-th/0303150 

[5] A. A. Belavin, A. M. Polyakov, A. B. Zamolodchikov, Infinite conformal symmetry in 
two-dimensional quantumfield theory, Nucl. Phys. B241 pp. 333-380 (1984) 

[6] L. F. Alday, D. Gaiotto, Y. Tachikawa, Liouville Correlation Functions from Four-dimensional Gauge 
Theories, Lett. Math. Phys. 91 pp. 167-197 (2010), 0906.3219 

[7] B. Eynard, N. Orantin, Algebraic methods in random matrices and enumerative geometry (2008), 
0811 . 3531 

[8] L. Chekhov, B. Eynard, O. Marchal, Topological expansion of the Bethe ansatz, and quantum 
algebraic geometry (2009), 0911.1664 

[9] L. Takhtajan, P. Zograf, Hyperbolic 2 spheres with conical singularities, accessory parameters and 
Kahler metrics on M(0,n) (2001), math/ 11217 

[10] I. Gelfand, L. Dikii, A Family of Hamiltonian structures related to nonlinear integrable differential 
equations (1978) 

[11] P. Bouwknegt, K. Schoutens, W symmetry, Adv.Ser.Math.Phys. 22 pp. 1-875 (1995) 

[12] A. Zamolodchikov, Higher equations of motion in Liouville field theory, Int. J. Mod. Phys. A19S2 pp. 
510-523 (2004), hep-th/0312279 

[13] H. Ooguri, A. Strominger, C. Vafa, Black hole attractors and the topological string, Phys. Rev. D70 p. 
106007 (2004), hep-th/ 4 0514 6 

[14] V. A. Fateev, A. V. Litvinov, Correlation functions in conformal Toda field theory I, JHEP lip. 002 
(2007), arXiv =07 9.380 6 [hep-th] 

[15] R. Dijkgraaf, C. Vafa, Toda Theories, Matrix Models, Topological Strings, and N=2 Gauge Systems 
(2009), 0909.2453 

[16] L. Chekhov, Logarithmic potential fi-ensembles and Feynman graphs (2010), 1009.5940 

[17] J.-E. Bourgine, Large N limit of beta-ensembles and deformed Seiberg-Witten relations, JHEP 1208 p. 
046 (2012), 120 6. 1696 

[18] A. Mironov, A. Morozov, A. Popolitov, S. Shakirov, Resolvents and Seiberg-Witten representation for 
Gaussian beta-ensemble, Theor.Math.Phys. 171 pp. 505-522 (2012), 1103.5470 



-25- 



